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SUMMARY 


Gas  phase  absorption  of  solar  radiation  has  been  suggested  as  an  energy  source 
to  power  space  propulsion  systems.  The  absorbing  gas  could  be  an  expendable 
propellant  or  could  be  used  in  a  closed  cycle,  exchanging  heat  with  an 
expendable  propellant.  This  study  was  conducted  in  order  to  develop  basic 
data  and  understanding  useful  for  the  evaluation  of  such  a  system  with  the 
absorbing  gas  being  a  halogen  or  an  interhalogen. 

Three  specific  chemical  systems  were  selected  for  study.  These  are  Cl 2 . 

IC1,  and  Cl^  *  H^.  In  this  latter  system,  the  chemical  energy  released  by 
the  photoini t iated  formation  of  HC1  from  the  elements  is  augmented  by  the 
radiation  energy  absorbed  by  Cl?. 

Experiments  were  conducted  in  which  photolysis  reactors  were  filled  with  the 
gases  over  a  range  of  pressures  and  flashlamps  were  fired,  thereby  exposing 
the  gases  to  a  short  pulse  (5  us  -  10  vis)  of  intense  radiation  intended  to 
spectrally  simulate  the  exoatmospheric  solar  radiation.  Optical  monitors 
provided  time  dependent  data  on  the  lamp  output  and  on  the  in-band  absorption. 
A  fast  pressure  transducer  provided  time  dependent  data  on  the  pressure  rise 
within  the  closed  reactors,  thereby  yielding  information  on  the  temperature 
rise  and  energy  absorption. 

Two  reactor  and  flashlamp  configurations  were  used.  The  first  was  a  cylin¬ 
drical  reactor  with  a  100  J  linear  flashlamp  mounted  on  its  axis.  This  setup 
provided  data  useful  in  demonstrating  energy  deposition  and  release  and  in 
validating  a  radiation-kinetics-thermodynamics  model.  Observed  pressure  rise 

was  limited  because  of  the  need  to  include  SF,  gas  in  the  reactor  to 

b 

alleviate  an  arcing  problem  and  because  of  the  relatively  low  ultraviolet 
output  from  the  lamp.  The  second  configuration  employed  a  9  kJ  linear  flash- 
lamp  and  a  small  reactor  adjacent  to  it  which  received  a  portion  of  the  lamp 
radiation.  This  lamp  was  richer  in  ultraviolet  output  and  the  reactor  did  not 
need  diluent  or  insulating  gas,  thus  leading  to  higher  observed  pressure  rises. 

Pressure  rise  data  in  Cl  and  IC1  corresponded  to  temperature  increases  of 
over  500K.  When  Cl_  +  H  were  exposed  to  the  flashlamp  radiaton, 


temperature  rises  of  up  to  2600K  were  observed.  In  the  Cl2  +  H2  system, 
radiation  augmentation  was  estimated  to  be  up  to  10-15  of  the  chemical  energy. 

An  extensive  kinetic  rate  package  was  assembled  for  the  reactions  studied  and 

this  data,  together  with  spectroscopic  data  on  the  flashlamp  and  on  the 

absorbing  gases  was  used  as  Input  to  a  Rocketdyne  numerical  analysis  model. 

Correlation  of  data  obtained  with  Cl2  +  SFfi  In  the  Internal  flashlamp 

configuration  was  quite  good.  Reasonable  correlation  of  data  obtained  with 

Cl.  +  H-  *■  SF,  was  achieved  after  significant  modifications  of  the  rate 
2  2  b 

set  to  include  chain  termination  by  SF,  and  0_.  With  the  SF  and  some 

b  2  b 

0^  impurity,  the  formation  of  HC1  from  the  elements  only  proceeded  part  way 
to  completion,  about  50%  at  most,  and  this  termination  was  successfully 
modelled  numerically. 

Analysis  of  the  data  obtained  with  the  external  flashlamp  arrangement  was 
quite  satisfactory  using  basic  spectroscopic  and  thermochemical  principles. 

Thus,  data  were  collected  and  analytical  methods  were  developed  which  have 
applicability  to  the  evalualton  of  solar  powered  space  propulsion  systems 
deriving  all  or  part  of  their  energy  through  gas  phase  absorption  of  solar 
radiation. 
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INTRODUCTION 
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Space  propulsion  systems  currently  in  use  or  planned  typically  operate  by 
expelling  high  velocity  exhaust  to  provide  thrust.  The  kinetic  energy  of  the 
exhaust  derives  from  potential  energy  which  is  stored  in  the  payload  as 
pressurization  or  as  chemical  energy  in  energetic  propellants.  It  is 
appealing  to  consider  advanced  space  propulsion  systems  in  which  all  or  part 
of  the  kinetic  energy  of  the  exhaust  is  provided  through  the  absorption  of 
solar  radiation  by  the  propellants. 

Advanced  space  propulsion  systems  based  on  these  concepts  must  be  evaluated  in 
terms  of  their  performance  relative  to  conventional  systems.  Evaluations  of 
this  type  have  identified  long -duration ,  low-thrust  missions  such  as  station 
keeping  and  attitude  control  as  potentially  attractive  applications  for 
advanced  propulsion  systems  deriving  energy  from  radiation  absorption.  One 
approach  to  converting  radiation  to  kinetic  energy,  which  is  the  approach 
addressed  by  this  study,  is  to  use  propellants  which  undergo  a  chemical  change 
upon  irradiation.  The  radiant  energy  is  thus  deposited  into  the  propellant  as 
the  heat  of  formation  of  the  new  species.  This  energy  is  subsequently 
recovered  and  converted  to  thrust.  In  some  systems  the  absorbed  radiation  may 
be  the  total  energy  available  while  in  other  systems  it  may  augment  chemically 
stored  energy. 

figure  1  shows  a  very  simplified  schematic  diagram  for  a  propulsion  system 
based  on  these  ideas.  The  system  shown  includes  a  fluid  flow  loop  with 
components  for  absorbing  solar  radiation,  for  transfering  the  energy  to  an 
expendable  propellant,  for  reconditioning  the  working  fluid,  and  for  pumping 
it  around  the  loop.  The  sketch  implies  a  closed  loop  for  the  fluid  which 
interacts  with  the  solar  radiation  and  a  separate  flow  system  for  the  expen 
dable  propellant,  although  variations  on  this  scheme  could  also  be  considered, 
furthermore,  the  source  of  energy  for  the  working  fluid  conditioner  and  pump 
is  not  identified  but  a  number  of  concepts  could  be  considered,  including 
tapping  some  of  the  energy  out  of  the  energy  transfer  device.  It  should  be 
emphasized  that  the  system  shown  here  differs  generically  from  the  solar 
dynamic  cycles,  such  as  are  being  considered  for  Space  Station  power  (Ref.  1) 
or  from  solar  heated  thruster  concepts  (Ref.  2)  in  that  all  or  part  of  the 
energy  added  to  the  working  fluid  is  directly  absorbed  via  an  overlap  between 
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Solar  propulsion  concept 


the  solar  spectrum  and  the  absorption  spectrum  of  one  or  more  components  of 
the  working  fluid.  The  work  reported  here  was  conducted  in  order  to  assist  in 
the  evaluation  of  solar  powered  propulsion  systems  of  the  general  type  shown. 

BACKGROUND 

Solar  Radiation 

The  spectrum  of  the  solar  radiation  outside  the  earth's  atmosphere  closely 

resembles  a  5900K  blackbody  spectrum  (Ref. 3  pp  3-34  to  3-29).  The  actual 

solar  spectral  irradiance  outside  the  atmosphere  is  well  represented  by  a 

standard  curve  proposed  the  Thekakara  (Ref .4)  and  shown  in  Fig.  2.  The 

coordinates  for  this  curve  are  tabulated  in  Ref.  3  and  4  at  0.005  and  0.010 

um  intervals  and  these  have  been  used  for  all  solar  spectrum  calculations 

discussed  here.  The  spectral  irradiance  integrated  over  wavelength  yields  the 

2 

total  exoatmospheric  solar  irradiance  of  1353  W/m  . 

Halogens  as  Solar  Radiation  Absorbers 

When  gas  phase  or  volatile  compounds  are  surveyed  for  photochemical  activity 
induced  by  solar  radiation,  the  halogens  and  their  compounds  show  up  clearly 
as  candidates.  The  diatomic  halogens  C 1 2 ,  Br2>  and  have  relatively 
broad  absorptions  peaking  at  about  330,  410,  and  490  nm,  respectively  (Ref.  5) 
These  are  shown  in  Figure  3.  The  quantum  yield  for  production  of  atoms  from 
these  molecules  is  essentially  one  from  250  to  500  nm,  falling  to  about  0.5 
for  1?  from  500  to  600  nm  (Ref.  6  pp  184-192).  This  leads  to  the  potential 
for  these  molecules,  taken  in  combination,  to  convert  approximately  30%  of  the 
solar  radiant  energy  into  the  heat  of  formation  of  ground  state  or  excited 
atoms.  When  diatomic  interhalogens  (Figure  4)  are  considered,  the  absorption 
range  is  extended  to  beyond  600  nm  (Ref.  7)  and  the  solar  absorption 
efficiency  to  above  35%.  There  are  other  more  complex  halogen-containing 
molecules,  perf 1 uoroaky 1 iod ides  for  example,  whose  photochemistry  is  well 
known  from  flash  pumped  laser  research  (Ref.  8).  Molecules  in  this  category 
were  not  considered  as  candidate  absorbers  because  of  potential  difficulties 
associated  with  byproduct  formation. 

Table  1  lists  some  information  on  the  halogens  and  interhalogens  relevant  to 
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their  consideration  as  potential  solar  absorbers.  This  Includes  the  photo- 
dissociation  reaction,  the  wavelength  range  absorbed,  the  wavelength  of  the 
peak  absorption,  and  the  percentage  overlap  of  the  absorption  spectrum  with 
the  solar  radiation  spectrum.  In  addition,  the  energy  absorbed  per  mole  of 
reaction  at  the  peak  absorption  wavelength  Is  given.  This  quantity  Is  a 
measure  of  the  energy  absorbing  capacity  of  the  halogen,  on  a  molar  basis. 
Further  data,  relevant  to  solar  energy  deposition  consideration,  Is  given  in 
Table  2. 

Recovery  of  Deposited  Energy 

There  are  a  number  of  ways  of  recovering  the  energy  deposited  in  halogens  by 
the  photolysis.  Two  recovery  paths  were  explicitly  considered  In  this  study. 
Irradiation  produces  halogen  atoms  according  to  the  reactions  of  Table  1  and 
also  excited  molecules  when  the  photolysis  yield  Is  less  than  1.  The  atoms 
are  typically  hot.  For  example,  the  peak  of  the  Cl2  absorption  at  330  nm 
provides  86.64  kcal  per  g  mol  of  C1?  dissociated:  this  is  28.79  kcal  more 
than  the  bond  energy.  The  first  recovery  path  considered  involves  systems 
containing  only  halogens.  After  thermallzation  of  the  hot  atoms,  recombina¬ 
tion  is  Induced  through  three-body  or  wall -catalyzed  recombination  reactions. 
In  the  case  of  pure  halogens,  the  parent  molecules  will  be  formed.  In  the 
case  of  interhalogens,  products  will  be  pure  halogens  and  diatomic  Inter¬ 
halogens.  Formation  of  more  complex  Interhalogens  has  been  assumed  negligible 
In  this  study.  The  diatomic  products  will  then  thermalize.  Thus  the  total 
photodeposlted  energy  will  be  recovered  as  a  temperature  and  pressures  rise  in 
the  gas.  Radiation  losses  in  this  type  of  system  have  been  assumed  negligible 
but  can  be  treated  using  a  thermal  radiation  feature  of  the  numerical  model. 

The  second  recovery  path  Involves  the  addition  of  a  fuel,  with  hydrogen  added 
to  chlorine  being  the  example  considered  here.  The  final  temperature  and 
pressure  of  the  product  gas  in  this  case  Is  a  result  of  thermal izat Ion  of 
photo-deposited  energy  plus  the  chemical  energy  provided  by  the  formation  of 
HC1  from  the  elements. 

The  well-known  exothermic  reaction  between  elemental  hydrogen  and  chlorine 
yields  44.15  kcal/g  mol  of  C 1 2  upon  weak  photoinitiation.  This  energy  Is 
augmented  by  up  to  86.64  cal/g  mol  of  Cl,,  for  complete  photodissociation  to 


ra 


& 


hot  ground-state  atoms.  Thus,  this  approach  yields  up  to  130.78  kcal/g  mol  of 
Cl?  dissociated  in  comparison  with  the  86.64  kcal/g  mol  of  Cl2  which  would 
be  available  with  no  hydrogen  present  or  the  44.15  kcal/g  mol  which  is  avail¬ 
able  with  weak  photoinitiation.  This  leads  to  higher  temperature  product 
gases  at  the  cost  of  added  complexity  and  lower  loop  efficiency  for  subsequent 
dissociation  of  the  HC1  back  to  H2  and  Cl2>  The  situation  is  shown  graphi¬ 
cally  In  Figure  5.  The  volumetric  energy  release  is  plotted  vs  the  wavenumber 
(and  wavelength)  for  the  three  situations  just  discussed.  At  room  temperature, 
the  bond  energy  for  Cl2  corresponds  to  20000  cm  1  (500nm).  Below  this 
energy,  the  photoinitiated  H2  +  Cl2  reaction  does  not  occur.  When  the 
chain  reaction  is  weakly  initiated,  a  negligible  amount  of  radiation  energy  is 
deposited  and  the  resulting  energy  release  is  independent  of  the  wavenumber, 
as  shown  in  the  lower  curve.  The  energy  deposited  in  Cl2  depends  on  the 

wavenumber  as  shown  by  the  middle  curve,  which  is  for  completely  dissociated 

Cl2-  When  completely  dissociated  Cl2  is  reacted  with  H2,  the  resulting 

energy  release  is  as  shown  by  the  top  curve,  which  is  the  sum  of  the  other 

two.  The  curve  for  energy  release  from  an  interhalogen  would  be  similar  to 
the  Cl2  curve  shown. 

OBJtCI  IVES  AND  SCOPE 


The  objective  of  the  work  reported  here  was  to  develop  data  and  analysis  for 
use  in  evaluating  the  feasibility  of  space  propulsion  systems,  such  as  the  one 
shown  earlier  in  Figure  1,  in  which  solar  energy  is  absorbed  via  photochemical 
reactions  and  ultimately  converted  to  kinetic  energy  of  the  exhaust.  In 
particular,  data  and  analysis  were  desired  so  that  the  solar  energy  deposition 
and  subsequent  thermal i zation  could  be  characterized.  Development  and  evalua¬ 
tion  of  detailed  schemes  for  implementing  the  generic  concept  shown  in  Figure 
1  were  not  included. 

The  class  of  absorbers  included  in  the  present  study  included  only  halogens. 
Three  categories  of  working  fluids  were  considered:  (1)  pure  halogens, 

(2)  mixed  halogens  and  (3)  halogens  plus  fuels.  Of  course,  in  an  actual 
system,  inert  diluents  might  be  added  to  any  of  these,  lhree  particular 
halogen -containing  systems  were  identified  for  study:  (1)  chlorine  gas, 

(2)  iodine  monochloride  vapor,  and  (3)  chlorine  gas  plus  hydrogen  gas.  In  the 
first  two  systems,  the  energy  released  originates  only  from  the  solar  energy 
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absorbed.  In  the  third  system,  the  energy  released  originates  from  the  solar 
energy  absorbed  plus  the  energy  associated  with  the  heat  of  formation  of  HC 1 . 
This  latter  chemical  system  had  been  discussed  previously  in  the  context  of 
its  utility  for  solar  conversion  (Refs.  9-13),  but  we  were  unable  to  find 
previous  reference  to  the  first  two  chemical  systems. 

APPROACH 

The  technical  approach  which  was  developed  for  meeting  the  objectives  included 
experiments  and  analysis.  Single  shot  flash  photolysis  experiments  were 
selected  because  of  their  relative  simplicity  and  the  ability  to  address  the 
full  volume  of  a  reactor  which  is  large  enough  to  make  wall  effects  relatively 
unimportant.  Also  important  was  the  consideration  that  xenon  flashlamps, 
which  were  selected,  provide  a  good  simulation  of  the  solar  spectrum,  espe 
daily  when  the  lamps  are  operated  at  high  power.  Figure  6  shows  that  the 
comparison  for  a  lamp  fired  at  5kJ,  0.75msec  is  quite  good. 

A  variety  of  methods  was  used  in  the  interpretation  and  analysis  of 
experimental  results.  However,  the  most  detailed  understanding  of  the 
chemistry  and  physics  of  the  energy  deposition  and  subsequent  reactions  was 
provided  by  a  Rocketdyne  model  designated  as  Radiation  Kinetics  Mixing  Model 
( RKM) .  The  model  was  used  in  the  time  domain  to  simulate  the  flash  photolysis 
experiment.  The  quality  of  the  correlation  between  the  data  and  computations 
was  quite  good  in  most  cases  and  validated  the  detailed  understanding  of  the 
basic  phenomena.  The  model  can  be  operated  in  the  space  domain  to  predict 
results  from  a  flowing  solar  reactor  but  this  was  not  done  in  this  study. 
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ison  of  solar  spectrum  and 
flashlamp  spectrum. 


EXPERIMENTAL  DESCRIPTION 


Ihe  experiments  which  were  conducted  involved  tilling  a  cylindrical  cell  with 
various  mixtures  of  He,  SFfc,  IC1,  Cl2  and  and  then  exposing  the 
contents  to  the  light  pulse  from  a  xenon  flashlamp.  The  principal  independent 
variables  were  geometry,  composition  of  the  gas  fill,  and  flashlamp 
intensity.  The  principal  dependent  variables  (observables)  were  all  recorded 
as  functions  of  time.  They  were  lamp  discharge  voltage,  lamp  discharge 
current,  transmitted  photon  flux  in  two  wavelength  bands,  and  pressure  rise  in 
the  cell.  The  experimental  arrangement  is  shown  schematically  in  Figure  7. 

The  experiments  were  conducted  at  the  Rocketdyne  Santa  Susana  Field  Labora¬ 
tories  (SSFL)  in  the  Chemical  Physics  Laboratory  (Building  315).  Three 
outdoor  test  cells  and  a  control  station  inside  the  blockhouse  were  used,  as 
shown  in  Figure  8.  Cell  3  contained  all  of  the  gas  cylinders.  The  rest  of 
the  gas  handling  system,  the  flashlamp,  and  the  photolysis  cell  were  in  cell 
2.  The  cold  traps  and  vacuum  pump  were  in  cell  1.  All  of  the  remote  control 
equipment  and  the  data  acquisition  were  in  the  blockhouse. 

The  apparatus  which  was  used  to  accomplish  each  of  the  functions  illustrated 
in  Figure  7  typically  evolved  during  the  program  in  response  to  perceived 
deficiencies  or  ideas  for  improvement.  In  the  following  paragraphs,  the 
apparatus  is  described  according  to  the  functions  shown  in  the  boxes  of 
Figure  7.  Only  configurations  which  produced  useful  data,  presented  later, 
are  described. 

GAS  HANDLING  SYSTEM 

The  gas  handling  system  was  a  vacuum  line  with  monel,  stainless  steel,  glass, 
and  Teflon  components.  The  system  is  shown  in  detail  in  the  flow  schematic  of 
Figure  9.  Components  are  identified  in  Table  3.  Operation  of  the  system  was 
quite  straightforward.  After  evacuation  of  the  system,  the  run  tanks  (11  and 
T2  in  Figure  9)  were  filled  to  the  appropriate  pressures  of  C 1 ^ ,  10%  C 1 ? 
in  He,  H_,  10%  H„  in  He,  He,  or  SF, .  H„  and  Cl0  were  never  mixed  in 
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this  operation.  This  filling  was  accomplished  with  one  operator  at  the 
cylinders  in  cell  3  (Ref.  Figures  8,  9)  and  one  operator  at  the  Switch  Pane1 
in  the  blockhouse  (Ref.  Figure  8).  The  operators  were  in  headset  communi - 
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TABLE  3. 

GAS  HANDLING  SYSTEM  COMPONENT  LIST 


Call 


Out 

(Ref. 

Mg.  9) 

Name 

GAS  B0T1LES 

B1 

Chlorine  in  helium 

bottle 

B1  .A 

Chlorine  bottle 

B2 

Helium  bottle 

B3 

Hydrogen  in  helium 

bottle 

83. A 

Hydrogen  bottle 

B4 

Sulfur  hexafluoride 

Bottle 

Rf.  GUI  61  OKI. 


R 1 

Cnlorine  in  Helium  regulator 

R1  .A 

Chlorine  regulator 

R? 

Helium  regulator 

R3 

Hydrogen  In  Helium  regulator 

R3.A 

Hydrogen  regulator 

R4 

Sulfur  hexafluoride  regulator 

VALVES 

VI 

Chlorine  in  Helium 
bottle  valve 

VI  .A 

Chlorine  bottle  valve 

Description/Type 


B-size,  10%  Cl 2  In  He,  maximum 
pressure  =  718  psi 

Size  3,  15  lb  liquid  C I 2 .  vapor 
pressure  85.3  psig  @  70  F 

K-size,  maximum  pressure  =  2400  psi 

B-size,  10%  Hp  in  He,  maximum 
pressure  -  700  psi 

Size  1A,  Gaseous  Hp,  800  psi 

Size  1A,  liquid  SF5,  vapor 
pressure  320  PSIG  @  70  F 

2 -stage  Matheson  H15C  Monel 
2 -stage  Matheson  H15C  Monel 
2 -stage  Matheson  3104 
2 -stage  Matheson  3104 
2-stage  Matheson  3104 
2-stage  Matheson  9-950 

Hand  operated 
Hand  operated 


VALVES  (continued) 


V2 

Helium  bottle  valve 

Hand  operated 

V3 

Hydrogen  In  Helium  bottle  valve 

Hand  operated 

V3.A 

Hydrogen  bottle  valve 

Hand  operated 

V4 

Chlorine  in  Helium  shut 
off  valve 

Hand  operated.  Monel 

V4.A 

Chlorine  shut  off  valve 

Hand  operated.  Monel 

VS 

Helium  shut  off  valve 

Hand  operated,  stainless 

V6 

Hydrogen  in  Helium  shut 
off  valve 

Hand  operated,  stainless 

V6.A 

Hydrogen  shut  off  valve 

Hand  operated,  stainless 

Cl 

Chlorine  in  Helium  check 
valve 

All  stainless 

Cl  .A 

Chlorine  Check  Valve 

All  stainless 

C-2 

Helium  Check  valve 

All  stainless 

C  -3 

Hydrogen  in  Helium  check  valve 

All  stainless 

C-3.A 

Hydrogen  Check  Valve 

All  stainless 

C4 

Surflur  hexafluoride 
check  valve 

All  Stainless 

V7 

Chlorine  purge 

Hand  operated,  stainless 

V8 

Hydrogen  purge 

Hand  operated,  stainless 

V9 

Chlorine  fill  valve 
(normal ly  closed) 

Air  operated,  stainless/Tef Ion 

VIO 

Hydrogen  fill  valve 
(normal ly  closed) 

Air  operated,  stainless/Tef Ion 

VI 1 

Chlorine  run  valve 
(normally  closed) 

Air  operated,  stainless/Tef Ion 

VI 2 

Hydrogen  run  valve 
(normally  closed) 

Air  operated,  stainless/Tef Ion 

VI 3 

Chlorine  Throttle  valve 

Hand  operated,  stainless  (9  turns) 

V14 

Hydrogen  throttle  valve 

Hand  operated,  stainless  (9  turns) 

u 


VALVES 

(continued) 

VI 5 

Bomb  valve 
(normally  closed) 

Air  operated,  stainless/Tef Ion 

VI 6 

Vacuum  valve 
(normally  closed) 

Air  operated,  stainless/Tef Ion 

VI  7 

Pump  valve 

Hand  operated,  stainless 

VI 8 

Bleed  valve 

Hand  operated,  stainless 

VI 9 

Chlorine  regulator 
purge  valve 

Hand  operated,  monel 

V20 

Chlorine  regulator  vent 
valve 

Hand  operated,  monel 

V21 

Chlorine  isolation  valve 

Hand  operated,  monel 

V22 

IC1  Fill  valve 
(normally  closed) 

Air  operated,  stainless/lef Ion 

V23 

Transducer  inlet  valve 

Hand  operated,  stainless 

V24 

Transducer  bypass  valve 

Hand  operated,  stainless 

V2S 

IC1  Shut-off  valve 

Hand  operated,  glass/Teflon 

V26 

Reference  vacuum  valve 

Hand  operated,  stainless 

V2  / 

Pump  bleed  valve 

Hand  operated,  stainless 

V28 

Sulfur  hexafluoride  bottle 
valve 

Hand  operated 

V29 

Sulfur  hexafluoride  shut 
off  valve 

Hand  operated 

INSTRUMENTS 

PI 

Chlorine  fill  pressure 

Heise  0-30  psig,  read  through 
blockhouse  window 

P2 

Hydrogen  fill  pressure 

Heise  0-60  psig,  read  through 
blockhouse  window 

I'.l 

Bomb  load  pressure 

Heise  0  60  psia,  read  through 
blockhouse  window 

T 


INSTRUMENTS  (continued) 


& 


Reference  pump  vacuum  gauge 


System  vacuir  jo  jge 


1C1  fITT  pressure 


Mixer  temperature 


Photolysis  reactor  temperature 


Transducer,  measure  cell 
pressure  Increase  after 
flashlamp  discharge 


Hastings  thermocouple  gauge, 
read  out  unit  at  pump 

Hastings  thermocouple  gauge, 
readout  unit  in  blockhouse 

Validyne  reluctance  manomater 
Model  DPI 5-43,  readout  unit  viewed 
through  blockhouse  window 

C.A.  thermocouple,  readout  in 
blockhouse 

C.A.  Thermocouple  readout  in 
blockhouse 

Model  S-102A15  Quartz  Transducer 
PCB  Piezotronics 


m 


cation  with  the  blockhouse  operator  able  to  view  gauges  indicating  pressures 
PI.  P2,  P3,  and  P5. 

The  gases  In  the  run  tanks  were  then  admitted  to  the  photolysis  cell  by  remote 

operation  of  valves  V-ll,  V-12,  V-15,  and  V-16.  Control  of  the  flow  to  the 

cell  was  accomplished  by  presetting  the  throttling  valves  V13  and  V14.  Back¬ 
ground  pressures  of  less  than  50  m  torr  were  achieved  in  all  parts  of  the 
systems  except  the  Cl2  run  tank  T1 ,  which  could  be  pumped  to  only  about 

300  m  torr.  The  pressure  gauges  PI,  P2,  and  P3  were  Bourdon-tube  type  gauges 

2 

manufactured  by  Helse.  The  had  either  0.1  or  0.2  lb/in  scale  divisions. 

The  gauge  used  to  read  P3,  the  photolysis  cell,  was  of  the  absolute  pressure 

2  2 
type  and  had  0.2  Ib/in  divisions  but  could  easily  read  to  0.1  lb/in  . 

The  zero  was  set  daily  and  an  error  of  <  1.5%  (actual  -  reading)  was  typically 

2 

observed  at  the  local  atmospheric  pressure  of  about  13.7  lb/in  . 

There  was  no  attempt  made  to  correct  the  pressure  readings  in  presenting  or 
Interpreting  the  data.  Based  on  this,  it  was  possible  to  add  gas  to  the  cell 
with  uncertainties  on  the  initial  and  final  pressure  readings  of  <  5  torr  and 
an  absolute  error  on  the  initial  and  final  pressure  of  about  1.5%. 

PHOTOLYSIS  CELL 

Two  different  photolysis  cells  were  used  during  the  program,  in  connection 
with  the  two  flashlamp  systems  which  were  used.  In  the  first,  the  flashlamp 
was  internal  to  the  cell  while  the  second  was  placed  adjacent  to  the 
flashlamp.  Both  cells  were  cylindrical  in  shape  and  of  stainless  steel 
construction. 

Internal  Lamp  Configuration 

The  final  configuration  of  this  cell  is  shown  in  Figure  10.  The  body  of  the 
cell  and  the  endcap  on  the  ground  end  were  of  stainless  steel  construction. 

The  endcap  on  the  high  voltage  end  was  made  of  Teflon  and  was  grooved  on  its 
inner  surface  to  alleviate  a  discharge  problem  noted  along  a  smooth  Teflon 
surface.  Vlton  0  rings  were  used  on  sealing  surfaces.  The  lamp  electrodes 
were  pressed  into  stainless  steel  clips  which  were  threaded  onto  stainless 
steel  rods  for  electrical  connection.  One  port  was  provided  for  evacuation 
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and  fill  as  shown.  Two  other  ports  were  used  for  optical  access  and  a  fourth 

port  was  used  for  the  pressure  transducer.  The  total  volume  of  the  cylindrical 

3 

portion  of  this  cell  was  731  cm  .  Other  volume  factors  included  the  flash- 

3 

lamp  at  about  5.5  cm  ,  and  the  ports  and  fill  system  volume  between  the  cell 

3 

and  last  valve  seat,  estimated  at  20-25  cm  .  The  dimensions  of  the  lamp 
capillary  were  12.7  cm  long,  0.3  cm  I.D.,  0.7  cm,  O.D.  Thus,  the  cylindrical 
area  over  which  light  was  admitted  to  the  cell  volume  at  the  flashlamp  outer 
surface  was  4.9  cm*\ 

External  Lamp  Configuration 


The  final  configuration  of  this  cell  is  shown  in  Figure  11.  The  cell  was  con¬ 
structed  of  heavy  wall  stainless  steel  with  Suprasil  windows  cemented  on  the 
ends  with  silicone  rubber  cement  (Red  RTV).  There  were  two  holes  through  the 
cell  sidewall  as  shown;  one  for  an  insulated  pressure  transducer,  and  one  for 


connection  to  the  gas  handling  system.  The  volume  of  this  cell  was  64.1  cm  . 


The  volume  of  the  delivery  line  from  the  cell  to  the  seat  of  the  last  valve 


was  estimated  to  be  less  than  3  cm  .  The  area  of  the  window  was  6.4  cm 


CELL  AND  LAMP  DIAGNOSTICS  AND  DATA  SYSTEM 


In  order  to  characterize  the  deposition  of  flashlamp  energy  into  the  gas  in 
the  photolysis  cell  and  the  subsequent  thermalization  of  that  energy,  certain 
time  dependent  measurements  were  made.  These  are  discussed  in  the  following 
paragraphs . 


Voltage  and  Current  Measurements 


The  voltage  waveform  across  the  small  flashlamp  was  measured  by  measuring  the 
voltage  across  a  2000:1  high  voltage  probe  which  was  connected  in  parallel 
with  the  lamp.  The  probe  was  a  LANL  design  rated  to  TOOkV.  Although  the 
2000:1  ratio  was  verified  using  lower  DC  voltages,  it  was  not  checked  up  to 
lOkV  on  this  program.  Therefore,  it  is  possible  that  discrepancies  noted 
between  observed  voltage  peaks  and  capacitor  charging  voltages  at  10  kv  could 
be  attributed  to  this  probe.  These  discrepancies  are  discussed  later  under 
small  flashlamp  system  electrical  characterization. 


24 


Photolysis  cell  for  use  with  external  flashlamp. 


Current  through  the  lamp  was  measured  by  observing  the  waveform  output  from  a 
100  A/V  current  transformer  manufacturer  by  Pearson  Instruments.  The  trans 
former  was  mounted  around  the  lamp  high  voltage  electrode  on  the  outside  of 
the  photolysis  reactor.  A  10X  attenuator  was  used  on  the  signal  Drior  to 
input  to  the  data  system. 

The  voltage  waveform  across  the  large  flashlamp  was  measured  by  recording  the 
output  from  a  895:1  voltage  probe  provided  by  the  vendor.  Current  measurements 
were  not  recorded  with  this  system. 

Optical  Measurements 

Two  optical  channels  provided  time-dependent  waveforms  of  the  flashlamp  pulse 
in  two  wavelength  bands.  Each  channel  consisted  of  a  quartz  fiber  bundle  with 
one  end  collecting  flashlamp  radiation  after  passage  through  the  reactor 
volume.  The  fibers  transmitted  the  light  into  the  control  station  and  onto 
silicon  photodiode  detectors.  Interference  filters  at  the  detector  end  were 
selected  to  pass  light  near  the  maximum  of  the  halogen  absorption  (X  =334  nm 
for  Cl^  and  X  -  500  nm  for  IC1)  for  one  channel  and  to  pass  light  beyond 
the  absorption  (X  >  550  nm  for  Cl?  and  X  =  650  nm  for  1C1)  for  the  other. 
Neutral  density  filters  (N.D.  =  3.0),  were  installed  on  the  long  wave  pass 
channel  to  avoid  detector  saturation  when  the  X  >550  nm  filter  was  used. 

The  systems  are  shown  schematically  in  Figures  12  and  13.  Optical  filter 
characteristics  are  shown  in  Figures  14  -  17.  Spectral  response  of  the 
detectors,  designated  as  UV  enhanced,  is  shown  in  Figure  18.  The  detectors 
were  used  with  either  the  3  K  Q  or  10  K  Q  feedback  resistor.  It  was  veri¬ 
fied  using  pulsed  LED  inputs  that  the  detector  rise  and  fall  time  were  well 
under  lus. 

Pressure  Measurements 

Dynamic  pressure  rise  in  the  photolysis  cell  following  flashlamp  discharge  was 
monitored  using  a  Model  S  102A15  quartz  transducer  manufactured  by  PCB  Piezo 
tronics.  Ihe  unit  had  a  natural  frequency  of  425  kHz,  a  rise  time  of  lus, 
an  average  response  of  23.1  mV/psi  and  a  linearity  of  <  1*  of  full  scale 
(lOOpsi  full  scale).  Ihe  linearity,  based  on  the  calibration  data  provided, 
was  much  better  than  that.  Calibration  was  provided  at  2  psi  intervals  to  10 
psi  and  at  20  psi  intervals  to  100  psi. 
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Difficulties  were  encountered  with  electrical  and  thermal  interf erances .  The 
electrical  problems  were  alleviated  by  mounting  the  transducer  in  a  nylon 
bushing  Instead  of  in  the  steel  reactor,  as  shown  earlier  in  Figure  11.  The 
thermal  problem  was  alleviated  by  covering  the  stainless  steel  transducer 
diaphragm  with  Teflon  tape,  thereby  providing  thermal  insulation  which  held 
off  the  thermal  shock  until  after  the  critical  pressure  rise  data  was 
collected. 

Data  Acquisition  System 

Time  dependent  signals  were  recorded  using  a  Model  4094  four  channel  digital 
recording  oscilloscope  manufactured  by  Nicolet.  Data  rates  up  to  0.5y.s  per 
point  on  all  four  channels  were  available.  A  Model  XF-44  dual-drive  floppy 
diskette  recorder  was  available  for  recording  the  data  from  each  test. 

flashlaMp  DRIVER  AND  FLASHLAMP 

Two  different  xenon  flashlamp  systems  were  used  in  this  program.  The  first 
was  based  on  a  small  flashlamp  mounted  inside  a  cylindrical  cell,  as  shown 
earlier  in  Figure  10.  The  second  system  was  based  on  a  larger  flashlamp  which 
was  placed  adjacent  to  the  cell  shown  in  Figure  11. 

SmaJ 1  Flashlamp  System  (~100  J  Discharge) 

The  driver  for  this  system  was  a  Model  457  A  Micropulser  manufactured  by  the 
Xenon  Corporation.  This  unit  was  equipped  with  capacitors  of  0.5  uF  and  2 
uF  which  could  be  charged  up  to  lOkV,  thus  providing  up  to  100J  of  discharge 
energy. 

The  flashlamps  were  also  manufactured  by  Xenon  and  carried  the  designation 
S-1084.  Ihe  capillary  section  of  the  lamp  was  12.6  cm  long  x  7mm  O.D. 

(3mm  I.D.).  The  flashlamps  were  filled  with  xenon  to  a  pressure  of  300  torr. 
Ihe  system  was  modified  as  follows  to  accommodate  operational  and  performance 
requirements  of  the  experiments. 


1.  Remote  operation  of  the  high  voltage  on-off  function  was  added  to 
satisfy  some  procedural  requirements. 


2.  A  remote  trigger  switch  circuit  with  a  synchronizing  pulse  output 
was  added  to  reduce  time  delay  jitter,  eliminate  switch  bounce 
effects,  and  reduce  electrical  interference  carried  by  the  synch 
pulse. 

3.  A  capacitor  of  200  nF  was  installed  across  the  flashlamp  to  provide 
a  higher  voltage  peak  at  the  beginning  of  the  discharge  to  improve 
lamp  ignition.  A  500  Q  bleed  resistor  was  installed  across  this 
capacitor. 

4.  The  braided  cable  on  the  ground  return  from  the  lamp  was  replaced 
with  2"  wide  copper  sheet  in  an  attempt  to  narrow  the  discharge 
waveform. 

A  trigger  wire  was  installed  for  all  of  the  lamp  firings  made  with  the  100  J 
flashlamp  system.  Several  hundred  lamp  firings  were  recorded  with  thi 
system,  more  than  half  of  them  at  100  J.  Reproducibility  of  the  electrical 
and  optical  characteristics  was  quite  good,  with  no  lamp  failures  encountered. 
The  overall  experimental  arrangement  for  the  experiments  conducted  with  this 
system  is  shown  in  Figure  19. 

large  Flash! amp  System 

The  large  flashlamp  system  was  manufactured  to  Rocketdyne  specifications  by 
the  Pulsar  Products  Division  of  Physics  International  Corporation.  The 
flashlamps  were  designated  as  Model  L-4258  and  were  manufactured  by  ILC 
Technology  Inc.  These  lamps  were  filled  with  30  torr  of  xenon,  had  an  0.0.  of 
24mm,  an  I.D.  of  19mm,  and  an  electrode  spacing  of  100  cm.  The  overall 
experimental  arrangement  for  the  experiments  conducted  with  this  sytem  is 
shown  in  Figure  20. 
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Figure  20.  Overall  arrangement  for  9kJ  flashlamp  experiments, 
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0BSERVA1 IONAL  RESULTS 


This  section  of  the  report  is  a  presentation  and  discussion  of  the  experi¬ 
mental  results  which  were  obtained  during  the  program.  The  discussion 
Includes  estimates  and  comments  on  the  quality  of  the  raw  data  and  methods 
used  to  reduce  the  data  to  physical  observables.  The  discussion  begins  with  a 
characterization  of  the  operation  and  performance  of  the  flashlamps.  Then 
energy  deposition  and  release  results  are  presented. 

FLASHLAMP  CHARAC1ER1ZA1 ION 


Electrical 


tion  in  the  Small  FlashlamD  System 


All  of  the  small  flashlamp  system  testing  discussed  in  this  report  was  done 

with  a  2  uf  capacitor  installed  in  the  Xenon  Corporation  Model  457A 

Micropulser.  Front  panel  meter  voltage  settings  of  5000,  7000,  or  10000  V 

2 

were  used,  giving  nominal  discharge  energies  (1/2  CV  )  of  25J,  49J,  and 
100J.  Early  in  the  testing,  difficulties  were  encountered  in  obtaining 
reliable  lamp  firings  due  to  shorting  or  arcing  inside  of  the  photolysis 
reactor.  A  number  of  steps  were  taken  to  alleviate  this  problem,  including 
the  addition  of  SF,  to  the  gas  fill.  It  was  discovered  that  at  least  100 

O 

torr  of  SF,  was  needed  to  suppress  the  unwanted  arcing.  Therefore,  all 

D 

tests  reported  with  this  system  included  100  torr  or  more  of  SF  In  the 

6 

photolysis  reactor.  The  role  of  this  gas  as  a  heat  sink  and  kinetics 
participant  was  not  particularly  desirable  and  moderated  observed  pressure 
rises  through  both  kinetic  and  thermal  effects,  as  discussed  later.  However, 
its  role  as  an  arc  suppressor  was  very  beneficial. 

Figure  21  shows  the  recorded  current  and  voltage  waveforms  for  a  10000V  lamp 
discharge  with  a  16  psia  fill  of  SF.  in  the  cell.  Also  shown  is  the  integral 
of  V(t).I(t)dt  performed  by  the  Nicolet.  '  The  spike  on  the  front  of  the 


(a)  N01E  FOR  INI ERPRE I  AT  LON  OF  GRAPH  SYMBOLS 

IHe  Nicolet  Digital  Oscilloscope  and  the  associated  Hewlett  Packard  x  y 
Digital  Plotter  were  used  regularly  in  recording,  manipulating,  and  displaying 
the  experimental  data.  The  graphs  produced  by  this  data  system  are  used 
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Electrical  character i s t i cs  of  10  kV  lamp  discharge. 


voltage  waveform  Is  caused  by  the  peaking  capacitor  discussed  earlier.  This 
spike  contributes  very  little  energy  but  helps  to  initiate  lamp  firing. 

The  waveforms  for  5  kV  and  7  kV  runs  are  similar  but  distinctive,  as  may  be 
seen  from  the  examples  in  figures  22  and  23.  ihe  principal  effects  on  the  V 
and  1  waveforms  of  increasing  the  voltage  and  energy  of  the  discharge  are  to 
narrow  the  waveforms  and  to  move  the  peaks  to  earlier  times,  with  the  larger 
effect  being  seen  between  5  kV  and  7  kV. 

Several  aspects  of  the  electrical  energy  deposition  performance  are  plotted  in 

figure  24,  which  is  a  composite  of  data  from  35  runs  made  on  December  12, 

1984.  The  parameters  plotted  are  the  peak  values  of  lamp  current  and  voltage, 

the  value  of  the  leading  voltage  spike,  and  the  time  integral  of  the  voltage- 

current  product,  which  is  the  energy  deposited.  Also  plotted  for  comparison 

are  the  set  voltage  from  the  Micropulser  panel  meter  and  the  stored  energy, 

2 

1/2  CV  .  based  on  this  value.  Data  for  each  parameter  are  plotted  as  a 
maximum,  minimum,  and  average  from  the  runs  for  each  voltage  setting,  with  11 
runs  at  5  kV,  8  runs  and  /  kV,  and  16  runs  at  10  kV  being  included. 

(a)  (footnote  cont'd. ).. .frequently  in  this  report.  This  note  is  intended  to 
facilitate  reading  these  graphs.  Four  plots  can  be  made  per  graph  and  the 
four  boxes  on  the  right  of  the  graph  provide  relevant  information.  The  right 
hand  edge  of  each  box  shows  the  type  of  line  used  to  plot  the  parameter,  with 
scaling  information  given  by  the  lettering  inside  each  box  as  follows: 

V/D  is  the  volts/division  on  the  vertical  axis  (normally  a 

calibration  factor  is  then  required  to  convert  to  engineering 
units) 

Vy  is  the  voltage  value  assigned  to  the  central  tick  on  the  vertical 
axis 

I/D  is  the  time  in  seconds  per  division  along  the  horizontal  axis 

11-  is  the  absolute  time  value  assigned  to  the  left  hand  end  of  the 

time  axis 

has  no  meaning  in  a  plot  of  a  waveform  vs  time. 


40 


waveforms  for  5,  7  and  10  kV  discharges 


waveforms  for  5,  7  and  10  kV  discharges. 
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Figure  2k.  Lamp  electrical  characteristics  (data 
from  12/12/8*0, 
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The  groupings  for  each  parameter  at  each  voltage  are  very  satisfying.  A 
larger  spread  noted  for  V  at  10  kV  is  caused  by  the  fact  that  the  spike 

Jr  l  Ivt 

narrows  at  this  voltage  and  the  0.5  us  data  sampling  rate  does  not  always 
catch  the  peak. 

V  is  about  80%  of  V  at  5kV  and  about  102%  of  it  a  10  kV.  Although 

rtAI\  jt  I 

the  V  measured  across  the  lamp  is  expected  to  be  slightly  less  than  the 

r  t  Al\ 

capacitor  charging  voltage,  the  large  deviation  between  the  two  voltages 
represents  a  discripancy.  Both  the  front  panel  meter  and  the  LANL  high 
voltage  probe  are  inherently  more  accurate  than  the  observed  deviaitons.  It 
is  believed  that  the  V(t)  data  was  Influenced  by  residual  ground  loop  problems 
which  were  quite  serious  with  this  system  and  which  had  been  reduced  but  not 
completely  eliminated,  The  problems  would  be  expected  to  be  more  serious  at 
higher  discharge  energies 

It  should  be  noted  that  the  data  grouping  shown  in  Figure  24  is  for  conditions 
ranging  from  16  psia  of  SF,  to  4  psia  of  SF,  plus  2  psia  of  10%  Cl„  in 

o  b  c 

He.  For  a  given  voltage,  with  SF  in  the  fill,  the  waveforms  are  nearly 

6 

identical  for  all  gas  fills.  Figure  25  shows  current  and  voltage  waveforms  at 
10  kV  for  the  extreme  SF.  containing  has  fills  just  noted.  It  may  be  seen 

D 

that  the  reproducibility  is  excellent. 

Optical  Characteristics  of  the  Small  Flashlamp  System 

In  this  section,  the  spectral  and  temporal  characteristics  of  the  flashlamp 
light  output  are  discussed.  The  observations  are  based  on  time  dependent  data 
collected  using  the  two  optical  channels  discussed  earlier.  Only  the  334nm 
bandpass  filter  (Figure  14)  and  the  550  nm  long  wavepass  filter  (Figure  15) 
were  used  with  the  small  flashlamp.  Figure  26  shows  the  output  from  these  two 
channels  together  with  the  voltage  and  current  waveforms.  The  optical  channel 
signals  are  designated  4>  n(t)  and  g>  (t).  Measurement  of  these 
waveforms  shows  widths  (FWHM)  of  5.8ps  for  the  current,  6.6  us  for  the  UV 
radiation,  and  8. Ops  for  the  visible  and  near  IR  radiation.  The  time  depen 
dent  integrals  for  electrical  energy  deposited  into  the  lamp  and  UV  and 
visible  light  out  are  compared  in  Figure  27.  Both  the  UV  and  the  longer 
wavelength  light  continue  to  be  emitted  after  the  electrical  energy  deposition 
is  complete,  with  the  long  wavelength  light  persisting  the  longest.  The 


25.  Reproducibil  i ty  of  waveforms  for  different  gas 


scharge  with  no  absorber 


no  absorber 


behavior  observed  for  the  UV  light  was  In  contradiction  to  expectations 
expressed  by  lamp  vendors  and  their  brochures.  There  was  some  concern  that 
the  waveforms  were  being  filtered  by  the  photodiode  circuits.  Tests  using 
pulsed  LEDs  as  sources  showed  that  the  rise  and  fall  times  of  the  output 
signals  were  well  under  one  microsecond,  thereby  alleviating  the  concerns  and 
leading  to  the  conclusion  that  the  waveforms  accurately  represent  the  light 
output.  Reproducibility  of  the  light  output  in  the  two  channels  was  quite 
good,  as  >ay  be  seen  from  the  waveforms  shown  in  Figure  28. 


As  expected,  the  light  output  variation  with  electrical  energy  input  was  not 

the  same  for  the  UV  as  it  was  for  the  longer  wavelength  radiation.  This  is 

shown  graphically  in  Figur°  29.  Time  integrated  output  signals  are  plotted  vs 

the  voltage  setting  for  4>LWp  and  for  <t>33Q.  The  ratio  of  these  quantities  is 

also  plotted.  As  the  discharge  electrical  energy  increases  by  a  factor  of  4 
2 

(from  1/2  CV  -  25  to  100  J),  the  long  wavelength  light  output  increases  by 
a  factor  of  2.8  while  the  light  output  near  330nm  increases  by  a  factor  of 
about  7.5. 

E lectrical  Energy  Deposition  in  the  Large  Flashlamp  System 

All  tests  done  with  this  system  were  conducted  with  the  29uF  capacitor 

2 

charged  to  near  25kV  (1/2CV  ^  9kJ).  A  current  monitor  was  not  available 
so  that  the  actual  energy  of  the  discharge  could  not  be  measured.  The  voltage 
as  determined  from  the  signal  out  of  the  895:1  voltage  probe  provided  inside 
the  unit  was  in  excellent  agreement  with  the  charging  voltage.  The  waveform 
of  the  discharge  voltage  at  a  charging  voltage  of  25  kV  had  a  width  ( FWHM)  of 
9.5  us . 

Optical  Characteristics  of  the  Large  Flashlamp  System 


Figure  30  shows  the  waveforms  for  the  discharge  voltage;  the  long  wavelength 


light,  4>|Wp;  and  the  ultraviolet  light,  <t>33Q . 


The  width  of  4>33Q  is  16  us 


while  the  width  of  is  22.5  us.  Both  the  UV  and  the  long  wavelength  light 

LWP 

appear  to  rise  at  about  the  same  rate  but  the  longer  wave-  length  light  takes 
longer  to  decay.  The  large  flashlamp  system  provided  significantly  more  short 
wavelength  light,  as  illustrated  by  Table  4. 
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Table  4.  FLASHLAMP  SPECTRAL  COMPARISON 


Lamp,  Voltage 


^330/^LWP^ 


Small,  5  kV 

0.064 

Small,  7  kV 

0.044 

Small,  10  kV 

0.18 

Large,  25  kV 

1.13 

(a)  This  Is  a  ratio  of  detector  signals  under  the  same  conditions  of 
attenuation  and  filtering,  not  a  ratio  of  optical  fluxes. 


Behavior  of  the  lamp  over  time  had  a  variation  as  shown  in  Figure  31,  where 
the  observed  signal  <&LWp  is  plotted  against  the  number  of  shots.  The  raw 
observed  peak  signal  is  plotted  as  well  as  the  observed  peak  adjusted  by 
variation  in  the  discharge  energy.  The  adjusted  values  appear  more  sensi¬ 
tive  to  the  voltage-dependent  adjustment  parameter  than  the  raw  data  are  to 
the  actual  voltage.  This  suggests  the  possibility  of  some  sort  of  current 
saturation  which  results  in  the  peak  light  output  being  relatively  independent 
of  voltage  at  settings  near  25kV.  There  appeared  to  be  a  slight  degradion  in 
light  output  during  each  test  day  and  a  possible  degradion  over  the  duration 
of  the  data.  This  latter  degradation  could  have  been  associated  with  varia¬ 
tions  in  the  fiber  bundle  positioning  from  day-to-day. 

Variation  in  the  UV  light  output,  as  indicated  by  <&  ,  was  larger,  at 

v«5U 

least  within  a  given  test  day.  This  may  be  seen  in  Figure  32,  which  is  a  plot 

of  the  ratio  of  the  $33Q  to  $LWp  signals  obtained  with  only  helium  in  the 

reactor.  In  the  case  of  the  8/14  data,  the  first  shot  produced  a  low  $33q  signal 

which  was  restored  after  window  cleanup  and  alignment  check.  Degradation  of 

the  UV  signal  during  this  day  is  tentatively  assigned  to  window  contamination 

associated  with  exposure  of  the  system  to  Cl^.  The  lines  drawn  on  Figure  31 

were  used  to  determine  the  factor  by  which  the  measured  <PLHP  should  be 

multiplied  to  determine  the  appropriate  I  for  each  shot  which  was 

3  3U  0 

made  with  Cl  in  the  reactor. 


52 


Summary  of 


NUMBER  OF  SHOTS  ON  LAMP 


T 


WPW 


FLASHLAMP  ENERGY  ABSORPTION  BY  Cl2 

When  the  gas  fill  In  the  reactor  contains  C 1 2 ,  part  of  the  flashlamp  energy 
is  absorbed  and  some  of  the  C12  is  dissociated.  As  the  thermalization 
described  earlier  proceeds,  the  pressure  in  the  constant  volume  reactor 
rises.  In  this  section,  the  observations  characterizing  the  absorption  and 
the  pressure  rise  are  presented  and  discussed.  The  results  are  presented 
separately  for  the  small  and  large  flashlamp  experiments  because  of  the 
differences  in  geometry  and  other  experimental  conditions.  In  both  setups, 
the  absorption  was  observed  via  the  optical  signal  and  the  pressure 
rise  was  observed  using  the  quartz  transducer. 

Ultraviolet  Absorption.  Small  Flashlamp  System 

As  C 1 ^  is  added  to  the  cell,  the  in-band  flashlamp  intensity  reaching  the 
detector  decreases,  as  shown  in  Figure  33.  When  the  cylindrical  geometry  of 
the  internally  mounted  flashlamp  is  analyzed  according  to  Beer's  law,  it  is 
found  that  the  ratio  of  the  observed  intensity  with  and  without  C 1 2  present 
should  be  given  by  the  familiar  exponential: 

110  _  in-epH 

I0(r)  -10 

where 

I ( r)  is  the  observed  Intensity  at  r  with  an  absorber  present 

IQ(r)  is  the  observed  Intensity  at  r  with  no  absorber  present, 

r  is  the  radial  position  of  the  measurement, 

e  Is  the  extinction  coefficient  in  1/mol. cm, 

p  is  the  Cl?  concentration  in  mol/1,  and 

1  is  the  radial  distance  from  the  lamp  surface  to  the 

measurement  position  in  cm. 

When  the  data  shown  in  Figure  33  are  compared  with  the  predictions  of  Beer's 
law,  the  results  are  as  shown  In  Figure  34.  When  the  actual  dimensions  are 
used  In  Beer's  law,  the  agreement  is  only  good  to  about  20-30%  on  Cl,,  pressure. 
If  the  absorption  data  are  used  to  derive  an  average  "effective"  path  length, 
the  agreement  is  quite  good  when  this  effective  path  length  is  then  used  in 
Beer's  law.  This  result,  with  the  effective  path  length  being  smaller  than 
the  actual  path  length,  is  consistent  with  light  reaching  the  detector  indi¬ 
rectly,  after  having  been  reflected  or  scattered  from  the  cell  walls.  Thus, 
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33.  for  various  Cl,  partial  pressures. 


Beer's  law  using  the  effective  path  length  can  be  used  to  determine  the 
chlorine  remaining  In  the  cell  after  a  run. 


Ultraviolet  Absorption.  Large  Flashlamp  System 

For  this  setup,  absorption  at  334nm  Inside  the  cylindrical  reactor  mounted 
adjacent  to  the  flashlamp  would  be  expected  to  follow  Beer's  law  if  the 
Incoming  wavefront  were  collimated  and  uniform.  Since  a  collimating  lens  was 
not  used,  strict  adherence  to  Beer's  law  was  not  expected. 

Data  from  the  tests  made  with  this  system  are  plotted  in  Figure  35.  Discre¬ 
pancies  between  observations  and  Beer's  law  approach  a  factor  of  two  for  the  5 
and  10  torr  experiment  3:  i.e.,  if  observed  I/IQ  and  Beer's  law  is  used  to 
infer  a  Cl^  pressure,  it  will  be  close  to  double  the  Cl^  pressure  added  to 
the  reactor  according  to  the  P3  pressure  gauge.  The  origin  of  this  discre¬ 
pancy  has  not  been  identified  and  it  is  probably  a  combination  of  factors, 
including  gauge  reading  uncertainty  and  optical  system  geometry  effects. 

Pressure  Rise  in  Small  Flashlamp  System 

All  tests  conducted  with  this  system  included  SF  in  the  reactor  to  prevent 

D 

arcing.  Figure  36  shows  a  typical  signal  from  the  quartz  transducer  on  two 
time  scales.  The  unexpanded  pressure  trace  before  and  after  the  event  is  seen 
to  carry  a  small  oscillation  with  a  period  of  about  17  ms.  This  apparently 
originated  in  the  transducer/power  unit  system  and  was  considered  acceptable. 
After  the  firing,  the  pressure  is  seen  to  rise  sharply  and  then  decay  back  to 
zero  in  about  300  ms.  The  decay  is  a  compound  phenomenon  associated  with 
actual  pressure  decay  and  the  transducer  characteristic  decay  time.  Sharp 
oscillations  with  a  period  of  about  5  ms  are  observed  during  the  decay.  When 
the  time  scale  is  expanded  64X,  as  in  the  second  trace  of  Figure  36,  oscilla¬ 
tions  with  a  period  of  about  1  ms  are  observed.  When  the  time  scale  is 
expanded  by  400X,  as  in  the  second  trace  of  Figure  37,  oscillations  with  a 
period  of  about  330  vs  are  observed.  All  of  these  oscillations  are  believed 
to  be  caused  by  acoustic  phenomena  associated  with  characteristic  dimensions 
of  the  experimental  apparatus. 


Figure  37  shows  the  optical  signal,  <^LWp  for  reference,  although  is  is 
largely  obscured  by  electrical  noise  on  the  pressure  trace.  This  noise  is 
damped  out  by  about  60  us.  Since  the  noise  has  been  damped  out  before  the 
first  pressure  wave  begins  to  arrive  at  the  transducer  at  about  250  us,  it 
was  considered  acceptable. 

A  series  of  tests  was  conducted  in  which  a  fixed  amount  of  SF,  and  a 

o 

variable  amount  of  C 1 2  was  added  to  the  reactor.  After  the  flashlamp 
discharge,  the  pressure  in  the  cell  was  observed  to  rise  rapidly  and  then 
decay.  The  observed  peak  pressure  values  are  shown  in  Table  5  The  peak 
pressures  were  taken  as  the  peaks  of  the  oscillations  shown  in  Figure  36  and 
37. 

The  results  from  a  series  of  tests  with  fixed  SF,  and  Cl_  and  variable 

O  C 

are  shown  in  Table  6.  Results  include  the  observed  peak  pressure  ise 
and  the  degree  of  Cl2  reaction,  which  was  determined  by  firing  the  lamp  a 
second  time  after  a  period  of  a  minute  or  less.  The  observed  ^33g  peak 
signal  was  then  used,  together  with  the  emperical  fit  shown  earlier  in  Figure 
33,  to  estimate  the  amount  of  of  Cl?  remaining.  As  seen  in  Table  5,  the 
Cl?  reacted  was  rather  low,  varying  from  near  10%  to  near  50%  for  the  condi¬ 
tions  shown.  This  is  discussed  further  in  the  analysis  section  of  the  report. 

Pressure  Rise  in  Large  Flashlamp  System 

Pressure  rise  signals  obtained  using  the  large  flashlamp  system  were  accom¬ 
panied  by  stronger  acoustic  phenomena  than  had  been  observed  with  the  100J 
system.  This  was  because  of  the  higher  energy  deposition  density,  smaller 
reactor  volume  and  location  of  the  transducer  with  its  diaphragm  nearly  flush 
with  the  reactor  wall.  Thermal  effects  were  also  observed  but  were  moderated 
by  covering  the  transducer  diaphragm  with  a  piece  of  adhesive-backed  TFE 
tape.  Figure  38  shows  the  pressure  trace  obtained  after  correction  for  a 
shifting  baseline,  with  620  torr  of  Cl2  in  the  reactor.  The  pressure  rises 
rapidly,  then  appears  to  follow  a  damped  oscillation  having  a  period  of  about 
900  us.  The  Fourier  transform  of  this  pressure,  shown  in  Figure  39,  contains 
not  only  a  fundamental  peak  near  1100  Hz,  but  also  several  weak  overtones. 

The  phenomena  leading  to  the  observed  pressure  rise  waveform  is  believed  to  be 
initially  a  travelling  wave  initiated  when  the  ultraviolet  light  is  deposited 
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Figure  38.  Pressure  trace  observed 


Into  a  thin  slab  of  C12  near  the  window.  This  probably  degrades  Into  more 
of  a  standing  wave  later. 

Table  7  displays  the  results  of  some  acoustic  calculations  and  shows  that  the 

observed  1100  Hz  Is  in  good  agreement  with  expected  acoustic  transit  times  for 

a  travelling  wave  at  300  K.  Closer  examination  of  the  waveform  on  an  expanded 

time  scale  (Figure  40)  shows  some  fine  structure  in  the  acoustic  pulses.  The 

first  few  oscillations  are  clearly  separated  into  two  peaks.  The  time  inter - 

-4 

val  between  peaks  1  and  3  of  7.8  x  10  s  should  correspond  to  one  round 

trip  of  a  distance  2  a  or  0.2  m.  This  leads  to  a  velocity  of  sound  of  256 

m/s.  Ignoring  y  variations,  a  temperature  of  430  K  can  be  derived  from 

this.  The  time  between  peaks  1  and  2  is  140  us.  The  ratio  of  this  time  to 
-4 

7.8  x  10  s  is  0.179,  in  good  agreement  with  the  ratio  of  the  distance 
between  the  transducer  and  the  near  window  to  the  total  reactor  length,  0.190 
(ref.  Fig.  11).  This  supports  the  idea  of  a  traveling  wave  bouncing  ,etween 
the  two  windows.  The  time  of  arrival  of  the  first  peak  after  the  light  pulse 
supports  a  speed  of  sound  near  280  m/s. 

A  number  of  approaches  was  tried  to  estimate  an  effective  pressure  and 
temperature  for  comparison  with  steady  state  analytical  results.  Figure  41 
shows  the  observed  A P  vs.  Cl?  pressure  for  the  AP  of  the  highest  peak 
and  for  the  average  AP  taken  between  the  peaks  of  the  first  and  second 
cycles  (e.g.  average  AP  between  peaks  1  and  3  of  Figure  40).  This  latter 
AP  is  believed  to  represent  an  effective  or  volumetrical ly  average  AP 
shortly  after  the  light  pulse.  Figure  41  also  shows  the  temperature  as 
derived  from  the  speed  of  sound  based  on  the  time  between  the  first  and  second 
peaks  (e.g.  peaks  1  and  3  in  Figure  40)  At  =  1/f  =  24/R.  Variations  in 
Y  were  ignored  in  deriving  the  temperature  from  c. 

When  and  C 1 ^  were  introduced  together  into  the  reactor  and  the  flash- 
lamp  fired,  the  pressure  rise  was  generally  much  higher  than  that  observed 
with  Cl,,  alone.  Figure  42  shows  a  typical  result  on  two  time  scales.  The 
fine  structure  observed  with  only  C 1 2  In  the  reactor  is  not  as  clear, 
probably  because  of  detonation  or  deflagration  phenomena  occurring  in  the 
reactor.  The  procedures  just  described  were  used  to  derive  AP  and  AT  from  the 
pressure  traces.  The  results  are  shown  graphically  in  Figures  43,  44,  and  45. 


pressure  trace  shown  in  Figure  38 
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Figure  41.  Observed  pressure  and  temperature  rise  in  C1_. 
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Figure  44.  Observed  pressure  rise  in  Cl^  ±  based 
on  average  of  first  two  peaks. 
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TABLE  7.  Acoustic  Frequencies  (Calculated) 


c=("yRT/M) 


f  =£ 

T  21 


2(a  +  8a/3tr) 


(a)  Expected  frequency  for  a  traveling  wave 

(b)  Resonant  frequency  for  a  cylindrical  colume,  where, 

c  =  speed  of  sound 
a  =  length 
a  =  radius 


FLASHLAMP  ENERGY  ABSORPTION  BY  IC1 

Whe  IC1  was  added  to  the  photolysis  reactor,  results  using  the  large  flashlamp 
system  were  obtained  which  were  quite  similar  to  those  obtained  with  Cl2. 

IC1  pressures  from  4.5  torr  to  29  torr  were  added  to  the  reactor.  Optical 
absorption  data  are  shown  in  Figure  46.  The  I/Iq  values  were  determined  as 
the  ratio  of  the  optical  signals  with  helium  and  with  IC1  in  the  reactor.  The 
500  nm  filter  (ref.  Figure  16)  was  used.  Figure  46  includes  the  results  of 
Beer's  law  calculations  for  wavelengths  corresponding  to  the  peak  filter 
transmission  (X  =  505  nm)  and  the  half  power  points  (X  =  488  nm  and  X  = 
5.24nm) . 

Pressure  rise  data  were  determined  as  described  earlier  and  are  plotted  in 
Figure  47.  Corresponding  Cl~  results  are  included  for  comparison. 
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Figure  46.  Observed  IC1  absorption  near  500  nm. 
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Mgure  47.  Observed  pressure  rise  In  IC1. 
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ANALYSIS  AND  INTERPRETATION 


This  section  of  the  report  Includes  results  of  analyses  conducted  to  correlate 
the  experimental  results  with  an  understanding  of  the  relevant  physics  and 
chemistry.  Analytical  methods  have  Included  a  variety  of  simple  calculations 
to  estimate  or  interpret  behavior  for  macroscopic  observables.  A  detailed 
kinetic  model  was  also  developed  and  the  relevent  coupled  equations  were 
solved  numerically  using  a  Rocketdyne  computer  program. 

NUMERICAL  ANALYSIS 

Computer  Program  Description 

The  computer  program  used  for  this  study  was  derived  from  the  Advanced  Cavity 
Code  (ACC),  which  had  been  developed  by  Rocketdyne  for  chemical  laser  cavity 
performance  evaluations  (Ref.  14,  15).  The  discussion  and  computations  pre¬ 
sented  here  apply  to  the  small  flashlamp,  cylindrical  geometry  arrangement. 
Modification  to  set  up  the  geometry  applicable  to  the  large  flashlamp  experi¬ 
ments  would  be  straightforward. 

Input  to  and  output  from  the  model  are  indicated  in  Table  8.  The  phenomena 
which  are  modeled  include  absorption  and  photodissociation,  two  body,  three 
body,  and  wall  reactions,  thermal  radiation,  heat  and  mass  transport  within 
the  gas  and  heat  transfer  with  the  walls.  Options  include  conducting  or 
insulating  walls,  and  the  addition  of  turbulence  to  enhance  transport 
phenomena,  and  specification  of  wall  reflectivity. 

For  photodissociation,  the  energy  transfer  is  considered  to  be  unidirectional, 
i.e.,  normal  to  the  flashlamp  surface  and  to  be  accompanied  by  photolytic 
reactions.  Photodissociation  is  treated  as  a  two-step  process.  The  first 
step  is  the  dissociation  of  molecules  into  atoms  in  which  the  photoenergy  is 
stored.  The  rate  of  dissociation  is  given  by  the  rate  of  photoenergy  absorbed 
divided  by  the  molecule  energy  for  dissociation  (i.e.,  photon  energy).  The 
second  step  is  the  deactivation  of  the  dissociated  atoms  to  ground  state. 

This  is  a  thermal izati or  process  by  which  the  stored  photoenergy  is  released 
and  heats  the  fluid.  The  local  photoenergy  flux  is  calculated  by  a  rate 


TABLE  8.  MODEL  OPERATION 


INPUT  TO  MODEL 

OUTPUT  FROM  MODEL 

LAMP  SPECTRUM,  /  0 (\,t)dt/At 

-[X^  (t,r) 

LAMP  TEMPORAL  BEHAVIOR,  /<*>(  \  ,t)  dX 

-P(t) 

-T(t,r) 

THERMAL  DATA 

"qwall 

KINETIC  RATE  PACKAGE 

SPECTROSCOPIC  DATA  PACKAGE 

REACTANT  ABSORPTION  COEFFICIENT 

PRODUCT  ABSORPTION  COEFFICIENT 

INITIAL  CONDITIONS:  P,  T,  [  X. ] 

GEOMETRY  AND  BOUNDARY  CONDITIONS 

Photo  flux 
Wavelength 
Time 

Concentration  of  species  i 

Pressure 

Temperature 

Heat  flux 


equation,  which  states  that  the  rate  of  decrease  in  energy  is  equal  to  the 
amount  absorbed. 


For  thermal  radiation,  the  energy  transfer  Is  random  in  direction  and  is  not 
directly  related  to  chemical  reacitons.  Thermal  radiation  in  the  gas,  due  to 
the  random  nature,  is  governed  by  Integro-differential  equations.  The  method 
of  differential  approximation  is  employed  by  which  the  equations  are  reduced 
to  ordinary  differential  form.  This  method  is  equivalent  to  the  Milne- 
Eddington  approximation  commonly  used  in  astrophysics.  The  radiative  heat 
flux  is  due  to  emission  minus  absorption.  The  emission  term  is  a  strong 
function  of  temperature,  while  the  absorption  term  is  itself  coupled  to  the 
radiative  heat  flux. 

Kinetic  Data  Set 


The  numerical  analysis  program  requires  as  input  a  chemical  kinetic  data 
package  including  species,  reactions,  and  rates.  The  reaction  rate  set  which 
was  developed  under  this  program  is  shown  in  Tables  9-11. 

Spectroscopic  and  Radiometric  Data  Set 

Spectroscopic  and  radiometric  data  input  to  the  RKM  model  comprises  parameters 
describing  flashlamp  intensity  as  a  function  of  wavelength  and  time,  reactant 
extinction  coefficients  as  a  function  of  wavelength  and  temperature,  and 
product  extinction  coefficients  as  a  function  of  wavelength  and  temperature. 

The  lamp  spectral  and  temporal  radiant  output  may  be  represented  by  <t>(\,t). 
Based  on  discussions  with  lamp  vendors,  the  total  radiant  output  was  taken  as 
35%  of  the  electrical  energy  input,  which  was  measured: 

4>(\,t)d\dt  -  0.35x  J  I(t).V(t)dt.  *  0.35  x  1/2  CV2 

where 

\  =  wavelength 
t  =  time 

I(t)  =  lamp  current 
V ( t )  -  lamp  voltage 
Vc  =  charging  voltage 
C  =  capacitance 


TABLE  11.  MOLECULAR  REACTIONS 


Cl  ♦  H2  ■  HCl(v)  ♦  H 
H  +  Cl2  ■  HCl(v)  ♦  Cl 


h'I(v)  ♦  KC1  *  PCI (v-1 )  *  MCI  1-6 
HCl(v)  +  H2  «  MCI (v-1)  +  H2  1-6 
HCl(v)  ♦  H  =  HCl(v')  ♦  H  V'< 

HCl(v)  ♦  Cl2  =  HCl(v-l)  ♦  Cl2  1-6 
HCl(v)  ♦  Cl  *  HCl(v’)  ♦  Cl 
HCl(v)  ♦  He  *  HCl(v-l)  *  He 
HCl(v)«HCl  M*HC1(v»1)vHC1(v'-1) 

Cl  +  IC1  «  Cl?  ♦  I 
Cl  -»  I2  -  1C1  ♦  1 
H  +  IC1  *  MCI (v)  ♦  I 


H  ♦  12  ■  HI  ♦  I 
Cl  *  HI  «  HCl(v)  *  I 


H  ♦ 

HI  • 

H2(v)  ♦ 

Cl 

*  C102 

*  C 1 2  ♦ 

H  ♦ 

C!0? 

o 

y~> 

X 

Cl 

>  H0? 

•  HC1(0) 

H  ♦ 

HOj  ♦ 

H2  4  o2 

-  8.250 

-  8.250 

-  8.250 

-  8.250 

-  8.250 

-  8.250 
>  -  .290 

-  7.569 
-14.126 
-44.561 


1.5  x  10'“ 
1.4  x  10'12 
1.4  x  ID'13 


2.3  x  10*14v 

5.3  x  10',5v 


7.5  x  10‘15v 
1.0  x  10-11 
4.9  x  10",7v 
2.2  x  10'J?(v*l) 


1.3  x  10  “v  2.6 

1.6  x  lO'20*  2.23 

8.3  x  10*11 

8.3  x  10  20v  2.0 

3.3  x  irii 

2.0  i  10‘29v  5.0 

6.6  x  10*J°  -1.0 


.3  x  10  20v I  2.0 


.0  x  10'29v  5.0 

.6,10^1  -’.0 


3.0  x  10" 
1.7  x  10* 

4.4  x  10* 

5.4  x  10* 


9.9  x  10* 


B.4  x  10" 


1.1  x  10'IU 
1.56  x  10-10 

r>  n  .  l n"  1 6 


Hew  Reaction  or  Hodified  Rate 


It  was  further  assumed  that  the  spectral  and  temporal  dependences  were 
uncoupled  and  were  given  by  normalizations  of  the  lamp  spectrum  F(X)  given 
by  Gonz  and  Newell  (Ref.  16)  and  shown  earlier  in  Figure  6  and  of  the  temporal 
behavior  observed  here  as  $33Q(t)  and  shown  earlier  in  e.g.  Figure  26. 

Thus,  <l>(X,t)  was  represented  as 


*(X,t) 


_  £1X) _ 

JF(X)dX 


^330^ 


x  0.35  x  J  l(t) 


or 


_  F(X) _ 

JF(X)dX 


*330^ 

•r<I>330^t^dt 


x  0 


.35  x  -2CV 


2 


x  V(t)dt 


Sample  Calculations 


The  computer  program  provides  graphical  output  in  the  form  of  parameter  values 

plotted  against  time.  Figure  48,  for  example  shows  the  pressure  In  units  of 
2 

dynes/cm  vs.  time.  Each  major  mark  on  the  time  base  is  1  ms  and  the  plot 

runs  to  10  ms.  Figure  49  shows  the  C 1 ^  concentration  in  mol  fraction 

plotted  against  the  same  time  base.  The  data  plotted  refer  to  node  15  in  the 

computation,  which  is  at  the  wall  of  the  cylindrical  reactor  with  the  small 

flashlamp  on  the  cylinder  axis.  The  computation  was  for  a  100  J  discharge  in 

207  torr  SF,,  62  torr  of  Cl.  and  248  torr  of  H_. 

b2  2 

Correlation  of  Results 


1 


Cl.  +  SF..  Results  from  a  series  of  tests  with  Cl.  and  SF.  in  the  small 
— 2 - 6  2  6 

flashlamp  system  were  shown  in  Table  5.  That  test  series  represented  a  con¬ 
stant  SF,  concentration  with  Cl.  concentration  being  varied  by  up  to  a 

V)  C 

factor  of  20.  When  the  experiments  were  modeled  using  the  rates  from  Tables 
9-11,  the  results  shown  in  Figure  50  were  obtained.  The  computed  peak  pres¬ 
sure  rise  increases  with  Cl?  concentration  in  excellent  agreement  with  the 
observed  trend.  A  relatively  constant  offset  of  about  2  torr  was  observed, 
with  the  computed  result  being  higher. 

Cl.  ±  H.  +  SF, .  Results  from  a  series  of  tests  with  Cl„,  H. ,  and 

£  c  b  2  2 

SF  in  the  small  flashlamp  system  were  shown  earlier  in  Table  6.  Obtaining 
o 

compu-  tational  results  in  agreement  with  the  Table  6  experimental  results 
proved  to  be  difficult.  When  the  difficulties  were  recognized,  it  was  decided 
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Figure  48.  Sample  computed  results»  pressure  vs  time. 


computed  results,  Cl  cc  .entration  vs  time 


Clo  PRESSURE  (torr  ) 


m 

n 


to  try  to  anchor  the  model  to  the  run  with  the  largest  observed  pressure  rise, 
3/11/85:  8/1.  The  difficulties  lay  In  getting  the  H2  +  Cl2  chain  reaction 

to  terminate  and  in  getting  simultaneous  agreement  with  the  observed  pressure 
rise  and  Cl2  consumption.  In  order  to  try  to  improve  the  agreement,  atom- 


atom  recombination  rates  with  SF,  as  the  third  body  were  varied.  Other 

o 

relevant  rates  were  considered  too  well  known  to  justify  changing  them.  In 
the  following  paragraphs,  attempts  to  Improve  agreement  between  experiment  and 
theory  are  summarized. 


Establishing  bounds  on  SF,  three  body  rate  excursions  was  accomplished 

0 

by  computing  3/11/85:  8/1  results  with  SF.  rates  from  1  to  100  times 

Q 

those  in  Table  9.  All  SF,  rates  were  multiplied  by  the  same  factor. 

0 


At  lOOx  the  baseline  rates,  much  less  reaction  occurred  than  had  been 


observed  experimentally.  At  lOx,  both  the  AP  and  the  AC1_  were  in 

-3  c  -3 

reasonable  agreement  with  experiment  in  the  3  x  10  s  to  6  x  1C  s 
time  regime.  However,  the  reaction  was  continuing  and  clearly  would  go 


to  near  completion  if  the  computation  were  continued.  This  revealed  an 
important  aspect  of  the  model  which  was  confirmed  in  subsequent  computa¬ 
tions:  SF,  alone  as  a  third  body  chain  terminator  could  not  account 
0 

for  the  experimentally  observed  AP  and  AC12  endpoints. 


2.  Estimating  the  best  choice  for  the  SF,  three  body  recombination  rates 

D 

was  accomplished  by  computing  3/11/85:  8/1  results  with  SF  rates  of 

O 

lOx,  30x,  40x,  and  50x  those  in  Table  7.  The  best  results  (AP  =  470 

_3 

Torr,  AC12  =  0.53,  at  t  =  6x10  s)  were  obtained  at  40x.  However, 
the  reaction  had  not  terminated  and  AP  and  AC12  were  still  rising. 

3.  Wall  reactions,  particularly  wall -catalyzed  recombination,  could  help  to 

terminate  the  chain.  They  were  included  at  this  point  (along  with  the 

40x  SF,  rates)  at  an  efficiency  of  100%.  A  few  percent  reduction  in 
o 

AP  and  AC1?  was  obtained  at  corresponding  times  but  the  chain  was 
not  terminated. 

4.  Some  turbulent  mixing  was  evidently  occurring  experimentally  as  evi¬ 
denced  by  acoustic  phenomena  In  the  cell.  The  laminar  diffusion  mixing 
was  enhanced  In  the  model  to  simulate  this  phenomena.  The  diffusion 
coefficient  was  Increased  by  about  lOOx  and  lOOOx  to  bring  the  reactants 
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to  the  catalytic  wall  more  rapidly.  Using  this  approach,  with  simul¬ 
taneous  variation  In  the  SFfe  rates,  good  agreement  with  AP  and 
ACl^  could  not  be  obtained.  Furthermore,  the  chain  was  not  terminated 
under  those  conditions  which  gave  the  best  AP,  ACl^  match  to  the 
experimental  data. 

5.  Two  other  phenomena  which  were  examined  were  the  H  +  SF  reaction  and 

o 

equilibrium  thermal  dissociation  of  Cl?.  Neither  of  these  effects 
could  account  for  the  observed  results. 

6.  Inclusion  of  0^,  resulting  from  air  Inleakage,  was  found  to  provide 

the  needed  chain  termination.  Figure  51  shows  the  sensitivity  of 

AC1  and  AP  to  0„  concentration  for  two  different  SF  rate 
2  2  o 

sets.  02  concentrations  in  a  reasonable  range  are  shown  to  provide 

agreement  with  the  observed  AP,  ACl^  results  for  the  test  case, 

when  coupled  with  SF  rates  10x-15x  these  In  Table  9.  More  impor- 

0 

tantly,  the  computations  showed  termination  of  the  chain  at  the  values 
plotted.  Wall  reactions  and  turbulent  mixing  were  not  Included  In  these 
runs. 

Figure  52  shows  the  results  of  the  model,  with  0.04  torr  of  02  and  the  lOx 
SF,  rate  set  applied  to  the  test  cases  of  Table  6.  The  pressure  rise  is 

D 

computed  to  be  somewhat  lower  than  the  experimental  values  except  for  the  test 
case  8/1.  Deviation  of  the  Cl 2  consumption  calculation  from  the  experiments 
is  quite  severe,  except  for  the  test  case.  At  this  point,  it  is  believed  that 
the  model  contains  all  the  features  needed  to  correlate  the  data  and  that 
optimization  of  the  agreement  with  the  data  points  shown  in  Figure  52  could  be 
accomplished  through  adjustment  of  the  SF^  rates  and  the  02  partial 
pressures  (which  were  not  necessarily  the  same  for  all  the  experiments). 
Because  of  the  lack  of  02  partial  pressure  data  and  the  expense  of  the 
computations,  this  optimization  was  not  attempted. 

ANALYSIS  OF  RESULTS  FROM  LARGE  FLASHLAHP  SYSTEM 

The  observational  results  presented  earlier  are  discussed  and  interpreted  In 
this  section.  Analytical  methods  involved  application  of  basic  spectroscopic. 
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chemical  and  thermodynamic  principles.  The  numerical  model  was  not  applied  to 
correlate  results  obtained  with  this  system. 


Energy  Deposition  and  Release  In  Cl 


2 


Flux  available  from  the  lamp  was  modelled  as  dependent  on  1/r  where  r  Is  the 

distance  from  the  lamp  centerline.  It  was  assumed  that  35%  of  the  electrical 

energy  deposited  was  avallabe  as  light  and  was  uniformly  distributed  over  the 

area  of  the  lamp  outer  wall  having  a  radius  of  1.2  cm  and  a  length  of  100  cm. 

The  flux  available  at  the  reactor  Inlet  window,  4.5  cm  from  the  lamp  center- 

2 

line,  was  thus  taken  as  1.1  J/cm  for  a  29pF,  25  kV  discharge  energy. 


The  energy  receiving  capacity  of  C12  In  the  reactor  depends  on  the  pressure- 
path  length  product.  If  a  334  nm  photon  (peak  of  Cl ^  absorption)  Is 
absorbed  per  molecule,  the  energy  capacity  In  the  10  cm  long  cell  depends  on 
the  C 1 p  pressure  according  to 

_ . . .  2.  .  P(torr) 

E(J/cm  )  =  57.4  x 

This  energy  capacity  of  the  Cl2  Is  plotted  In  Figure  53.  When  the  fraction 
of  the  lamp  energy  at  each  wavelength  (from  the  assumed  lamp  spectrum)  is 
multiplied  by  the  Cl2  absorption  as  a  function  of  pressure  and  the  product 
is  Integrated  over  the  wavelength  range  of  the  Cl2  absorption,  the  amount  of 
lamp  flux  which  Is  absorbed  Is  also  shown  In  Figure  53.  The  ratio  of  the 
energy  absorbed  and  the  energy  capacity  gives  the  fraction  of  the  Cl2 
absorbed,  which  Is  also  plotted  In  Figure  53. 


The  energy  absorbed  by  the  Cl2  leads,  after  thermal izat ion,  to  a  temperature 
rise.  This  can  be  estimated  using  the  amount  of  energy  absorbed  and  the  heat 
capacity,  according  to 


AT 


E 

nC 

v 


2 

where  n  =  number  of  mols/cm  and  C  is  the  heat  capacity  at  constant 

-l  v  -1 

volume,  which  Is  6.132  cal  mol  deg  at  300K. 


ENERGY  (J/cm2).  AT  (K),  DISSOCIATION  FRACTION 


Figure  53.  Analysis  for  temperature  rise  In  C1^ . 
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The  AT  estimate  Is  also  plotted  In  Figure  53.  For  comparison,  observational 
results  presented  earlier  are  also  plotted  on  the  figure.  The  AT  results 
shown  earlier  in  Figure  41  are  plotted  directly.  The  AP  results  shown  in 
Figure  41  were  used  to  derive  a  AT  using 

AP/P  =  AT/T 

with  T  taken  as  300K. 

Considering  the  assumptions  made  In  deriving  the  AT  predictions  in  this 
section,  and  the  difficulties  described  earlier  In  obtain  AP  and  AT  values 
from  the  data  records,  the  agreement  Is  quite  satisfactory.  The  fact  that  the 
different  AT  determination  methods  provide  agreement  with  the  prediction  in 
different  pressure  regions  may  provide  further  Insight  into  the  predictive  and 
data  reduction  methods  but  this  has  not  been  pursued. 

Energy  Deposition  and  Release  in  Cl^  » 

The  final  temperature  in  this  system  is  achieved  through  the  thermal izat ion  of 
the  absorbed  flashlamp  radiation  plus  the  chemical  energy  released  by  the  exo¬ 
thermic  formation  of  HC1  from  the  elements. 

-T  _  E(absorbed)  »  E(chemlcal) _ 

nHClCv,HCl  *  n<H2Cv,H2  +  n'ci2Cv,Cl2 

where  n'H  and  n'cl  represent  any  reactants  which  remain  after  the  chain 

reaction  has  been  terminated.  For  the  analysis  presented  here,  complete 
reaction  Is  assumed.  With  the  10.3  torr  and  20.7  torr  partial  pressures  of 
Cl2  for  which  Cl2  *■  H2  tests  were  conducted,  the  absorbed  energy, 
discussed  In  the  previous  section.  Is  about  10*  of  the  chemical  energy.  This 
percentage  goes  down  as  the  Cl2  pressure  goes  up,  resulting  in  a  slightly 
lower  predicted  AT  for  the  20.7  torr  than  for  the  10.3  torr  cases.  Figure 
54  shows  the  calculated  temperature  rise  plotted  as  a  function  of 
ratio.  Experimentally  determined  results  are  included  on  the  plot  for 
comparison.  The  experimental  values  are  those  determined  using  the  average 
pressure  rise  over  the  first  two  AP  peaks  and  then  employing  AP/P  =  AT/T 


O'. 


Figure 


to  determine  AT  values.  Agreement  between  the  calculated  and  experimental 

results  are  quite  satisfactory  at  high  H2 : C 1 2  rat*os-  As  rat^°  is 

reduced,  the  reaction  to  form  HC1  does  not  go  to  completion  and  the  chemical 

energy  Input  Is  reduced,  leading  to  a  lower  observed  AT.  The  temperature 

rise  calculations  shown  In  Figure  54  were  obtained  using  2  or  3  Iterations  on 

the  specific  heat.  Initial  AT  values  were  calculated  using  Cyat  300K. 

The  calculations  were  repeated  using  a  C  value  for  an  average  (T  *■  T00)/2 

V  r INAL 

temperature.  The  procedure  was  repeated  again  If  necessary. 


CONCLUSIONS 


The  work  reported  here  was  conducted  to  obtain  and  demonstrate  a  basic  under 
standing  of  the  phenomena  of  energy  deposition  and  release  in  halogens,  mixed 
halogens,  and  halogens  plus  fuel.  The  results  are  useful  for  evaluation  of 
space  propulsion  systems  which  derive  all  or  part  of  their  energy  from  solar 
absorption.  However,  the  scope  of  work  did  not  include  engineering  or  system 
evaluations  of  such  space  propulsion  systems.  Therefore,  the  conclusions 
presented  here  relate  to  basic  phenomena  rather  than  potential  feasiblity  of 
applications. 


K 


The  flash  photolysis  approach  to  evaluation  of  energy  deposition 
and  release  provides  an  easy  and  convenient  way  to  screen 
reactants  and  reactor  operating  conditions.  Conditions  which 
provide  high  energy  absorption  and  release  can  be  quickly  identi¬ 
fied.  Quantitative  determination  of  the  energy  absorption  and  of 
energy  release  contain  uncertainties  due  to  acoustic  phenomena 
and  incomplete  information  regarding  flashlamp  output. 

Experimental  determination  of  energy  deposition  and  release 
provided  self  consistent  data  which  was  in  satisfactory  agreement 
with  analytical  predictions  and  correlations. 

Simple  analysis  based  on  spectroscopic  and  thermochemical  princi¬ 
ples  can  provide  a  satisfactory  correlation  of  pressure  and 
temperature  increase  induced  by  absorption  of  radiation.  However, 
the  analysis  does  not  predict  the  degeee  to  which  the  chain 
reaction  goes  to  completion. 


The  Rocketdyne  numerical  analysis  program.  Radiation  Kinetics 
Mixing  (RKM)  provides  good  predictions  of  energy  deposition  and 
release  for  halogen-only  systems.  For  systems  involving  a  chain 
reaction,  the  role  of  impurities  and  diluents  as  chain  termi¬ 
nators  can  be  explored.  By  modification  of  the  chain  termination 
rates,  the  model  appears  to  be  capable  of  correctly  predicting 


96 


the  time  evolution  of  properties  in  the  photolysis  reaction.  The 
spatial  evolution  feature,  applicable  to  a  flowing  reactor,  was 
not  applied  here.  It  is  believed  that  the  computations  performed 
are  the  most  detailed  to  date. 

The  kinetic  rate  package,  which  was  developed  using  direct 
literature  values  and/or  analogies,  appears  to  be  generally 
satisfactory  for  the  chemical  reactions  studied.  Some 
uncertainties  remain  regarding  recombination  rates  involving 
SF,  as  a  third  body.  The  kinetic  data  set  is  believed  to  be 

D 

the  most  comprehensive  and  accurate  available  for  the  chemical 
rations  studied. 

Energy  deposition  in  Cl,,  and  in  IC1  initially  at  room 
temperature  can  lead  to  a  temperature  rise  of  up  to  nearly 
1000K.  This  corresponds  to  nearly  10*  of  the  halogen  being 
dissociated.  Enchanced  absorption  coefficients  at  higher 
temperature  would  increase  the  energy  depositions  but  this 
phenomenon  has  not  been  explored  in  this  program. 

Photoinitiated  reaction  of  Cl?  +•  H?  mixtures  can  lead  to 
temperature  rises  of  a  few  thousand  degrees  K.  As  much  as  10-15* 
of  the  energy  released  can  originate  in  absorption  of  radiation. 

The  photoinitiated  reaction  between  Cl?  and  H2  fails  to  go  to 

completion  at  low  pressure  without  excess  H^.  When  diluents  or 

Impurities  are  present,  these  terminate  the  chain.  SF,  is 

6 

moderately  effective  as  a  chain  terminated  and  0?  is  extremely 

effective.  In  the  reactions  observed  here  with  SF  diluent  and 

6 

0,,  Impurity,  the  highest  Cl?  utilization  was  only  about  50*. 


RECOMMENDATIONS 


Based  on  the  results  and  understanding  developed  In  this  program.  It  would 
appear  logical  and  beneficial  to  pursue  this  general  area  of  research  and 
development.  The  recommendations  presented  here  are  of  three  general  types. 
The  first  type  represents  engineering  analysis  In  order  to  assess  the  feasi¬ 
bility  of  space  based  solar  propulsion  systems  employing  the  principles 
discussed  here.  The  second  type  represents  additional  experimental  and 
analytical  research  leading  to  Improving  the  data  base  and  understanding 
relative  to  solar  converions  via  gas  phase  absorption  of  solar  radiation.  The 
third  type  of  recommendation  deals  with  development  of  proof  of  principle 
solar  reactors  applicable  to  space  propulsion. 

1.  Engineering  analysis 

It  is  recommended  that  a  small  design  analysis  study  be  conducted  to 
configure  and  optimize  space  propulsion  systems  based  on  halogen  gas  phase 
absorption  of  solar  radiation.  The  study  could  begin  Immediately  using 
the  results  and  understanding  developed  in  the  present  program  plus  other 
open  literature  data. 


2.  Additional  experiments  and  analysis. 

It  is  recommended  that  the  experimental  and  analytical  approaches 
developed  in  the  present  program  be  applied  to  evaluation  of  some 
Important  aspects  of  solar  radiation  by  halogens.  Enchancement  of  the 
absorption  through  heating  of  the  reactants  and  through  use  of  mixed 
halogens  should  be  explored.  Improvements  in  the  experimental  approach 
through  improved  characterization  of  the  flashlamp  output  and  measurement 
of  total  energy  absorption  should  be  included. 

3.  Small  scale  reactor  development. 

It  is  recommended  that  a  small  scale  solar  absorber  be  developed  and 
tested  using  terrestrial  solar  flux  .  The  device  should  be  of  a  flowing 
gas  closed  loop  design  to  allow  some  of  the  absorbed  energy  to  provide 
system  heating,  therby  enhancing  absorption  coefficients.  This  heating 
could  also  be  simulated  through  electrical  heating  of  a  flowing  open  loop 
device.  Instrumentation  would  provide  for  controlled  extraction  and 
measurement  of  the  deposited  energy. 
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